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ABSTRACT 

In the epoch of precise and accurate cosmology, cross-confirmation using a variety of cosmographic 
methods is paramount to circumvent systematic uncertainties. Owing to progenitor histories and 
explosion physics differing from those of Type la SNe (SNela), Type Il-plateau supernovae (SNell-P) 
are unlikely to be affected by evolution in the same way. Based on a new analysis of 17 SNe II-P, 
and on an improved methodology, we find that SNe II-P are good standardizable candles, almost 
comparable to SNela. We derive a tight Hubble diagram with a dispersion of 10% i n distance, using 
the simple correlation between luminosity and ph otospheric ve l ocity introduced by iHamuv fc Pinta 
12002). We show that the descendent method of Nugent et all ()2OO60 can be further simplified and 
that the correction for dust extinction has low statistical impact. We find that our SN sample favors, 
on average, a very steep dust law with total to selective extinction Ry < 2. Such an extinction 
law has been recently inferred for many SNela. Our results indicate that a distance measurement 
can be obtained with a single spectrum of a SNII-P during the plateau phase combined with sparse 
photometric measurements. 

Subject headings: cosmology: observations — distance scale — dust, extinction — supernovae: general 

1. INTRODUCTION 

Compelling evidence for cosmic accelerat i on comes from distance measurenients t o T ype la supernovae (SNela; e.g., 
lAstier et al.] 120061 IWood-Vasev et al.ll200ii : iRiess et all 120071 : iKowalski et"ani2008L see Il^ilippenkoll2005b l for a review 
of earlier studies). Yet even with unlimited observational resources to improve statistical uncertainty, SN la cosmology 
ultimately faces systematic uncertainties stemming from an incomplete physical model of the phenomena. Indeed, 
since the physical conditions which give rise to the progenitors of the events are extremely difficult to discern, with no 
real consensus having been reached thus far, determining the nature of the evolution of SNela at high redshift remains 
a crucial challenge for precise and accurate measurements of the fundamental cosmologi cal parameters. While several 
mechanisms have been devised to measure, test, and constrain SNIa systematics fe.g.. lEUis et all 120081 : [Foley et all 

■ 12008*), the possibility of source redshift evolution warrants a continued exploration of additional, complementary 
cosmographic methods. 

A simple observation provides a compelling foundation for the use of Type Il-platea u supernovae (SNe II-P) as 

■ stand ardizable candles: the progenitors of several SNe II-P have now been uncovered (e.g.. lLi et al.ll2067tlSmartt et al.l 
l2008l and references therein), and all are found to be red supergiants with somewhat similar initial stellar masses 
(typically 8-16 Mq). Consequently, the environments and pre-SN evolution of such objects are reasonably tractable. 
Having undergone core collapse and envelope ejection, these red supergiants become SNe ob servationally define d by 
the presence of hydrogen in the spectra and a "plateau" phase in their light curves (e.g.. iBarbon et aT]ll979l see 
lFilippenkolll997l for a review). SNe II-P have been detected to redshift z « 0.5, and any evolutionary trends are likely 
to differ from those of SNe la. 

As in the case of SNe la, the intrinsic inhomogeneity in SN II-P peak luminosity, and that produced by extrinsic 
factors such as dust, can be calibrated. There are two different approaches to determining distance: the theory- 
based method, and the empirical stand ardized-candle method. The theoretical appr oach branches into two: (1) the 
expanding photosphere method (EPM: iKirshner fc Kwanl ll975l : [Eastman et al.|[l996l ). a historical descendent of the 
Baade-Wesselink met hod for variable sta rs (Baade 1926), and (2 ) the more modern "synthetic spectral atmosphere 
fitting" method (e.g., iBaron et all 120041 : iDessart & Hillier |2006[ ). Both techniques require the comparison of high 
signal-to-noise ratio (S/N) photometry and spectra of a given object to model atmospheres; for the second method, 
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synthesized spectra are computed for each SN in detail. Distances based on improved versio ns of the theoretical 
modeling for a handful of well-observed SNell-P have been shown to be precise to 10% (e.g., [Leonard et al .1120031 : 
lHamuvl[2005l : [Baron et a l. 2004; Dcssart & Hillicr 20061: [Pessart et ani2008f ). not greatly inferior to the best precision 
achievable with SNela (7%;|Astier_et_alJ|2006). The requirement of high S/N observations, however, limits the efficacy 
of this approach for events at larger distances. 

The standardized candle method first suggested bv lHamuv fc Pintol (|2002l : hereafter HP02) provides an independent 
empirical way to achieve distances to SNell-P, and it is strongly anchored in simple physics. In more luminous SNe, 
the H-recombination front is maintained at higher velocities; the photosphere is farther out in radius. When a SNII-P 
is on the plateau phase of its light curve (which lasts for around 100 rest-frame days), there is a strong correlation 
between the velocity of the weak Fe II lines near 5000 A ~ that trace the photospheric velocity - and the luminosity on 
the platea u. Extinction c orrections are based on a variety of different methods. This empirical correlation was further 
studied bv iHamuvl (|2005l : hereafter H05); application to 24 SNell-P in the Hubble flow yields a Hubble diagram in 
t he / band with a sca tter of 15% in distance. 

iNugent et al.l (|2006t hereafter N06) modified this technique by simultaneously combining both the extinction correc- 
tion (using the rest-frame V — I color on the plateau) and the Fe II velocity correction to arrive at a simple correlation 
between these parameters and luminosity. Given the smaller quantity of the data required, the NOG method, or a 
variant thereof, is probably the only framework that could be followed cost-effectively at high redshifts. 

The number of SNell-P for which distances have been derived, using any available method, is ~20. The paucity of 
such objects, and the intractable biases that emerge from the way the sample was constructed, warrants an analysis 
of more SNe, and push for improved and better-tested correlations that will allow one to achieve the precision and 
control of systematics that could make these SNe competitive tools for precision cosmology. 

In this work, we add 17 new objects to the current sample of low-z SNell-P, almost doubling the size, and reanalyze 
3 previously published events. We discuss the full sample construction in Sj2l In S}3] we present and modify the N06 
fitting method. We discuss in blithe sample culling, derive a new Hubble diagram, and find that the scatter is '--'10% 
in distance; we also consider the robustness and possible biases in the method, and compare distances derived for SNe 
that occurred in the same host galaxy, checking for internal consistency. Our conclusions are summarized in fj5l 

2. THE SUPERNOVA SAMPLE 

2.1. The Previously Existing Samples 

NOG compiled a list of 24 SNe II-P, 6 of which had n ot yet been pubUshed ; one new event (SN 2004dh) was studied 
as part of the Caltech Core Collapse Project (CCCP: fCal-Yam et al.ll2007h and is at low re dshift. while the o ther 5 
are at moderate redshifts, 0.1 < z < 0.3, found in the course of the Supernova Legacy Survey (jAstier et al.ll200(j ). The 
other 18 SNe, previously collected by H05 and HP02, were found and observed between 1986 and 2000, with a variety 
of instruments. All of those SNe had z < 0.06, though the majority were at z < 0.01. For such nearby objects peculiar 
velocities of the host galaxies dominate the error budget, thus complicating the derivation of robust correlations and 
distances. 

Since a "standard" Hubble diagram of distance modulus vs. redshift assumes that the sole contribution to the 
redshifts is cosmological, we apply corrections to the heliocentric redshifts in order to account for the peculiar velocities 
of nearby galaxies. Following previous authors, in the analysis bel ow we correct the h eliocentric redshifts derived from 
the SNe and their host galaxies to the local Hubble fiow using the lTonrv et al.l ()2000[ ) prescription. We further assume 
a Hubble constant Ho = TOkms^^Mpc"^. Since we later fit for the absolute magnitude of the SNe, Hg merely 
serves as a convenien t scaling const a nt. F or SN 1999em, we derive a redshift using the Hubble law and the Cepheid 
distance measured bv [Leonard et al.l ()2003D . 11.7± 1 Mpc. We re-reduced the images of SN 2004dh and report here new 
photometry for that event. Finally, following H05, we remove SN 2000cb from the sample, since it is now believed to be 
a SN 1987A-like event rather than a classic SNII-P. Hereafter we will refer to these 23 SNe as the "HPN sample" (for 
Hamuy, Pinto, and Nugent). 

2.2. The KAIT Sample: Nearly Doubling the Low-Redshift Census 

The 0.76-m Katzman Automatic Imaging Telescope (KA IT) is m ostly known as the SN discovery engine of the 
Lick Observatory SN Search (LOSS; Filippcnko et al. 2003: [FTlippe nko 2005^, but a substantial fraction of its time 
is dedicated to follow-up broadband photometry of nearby SNe. As part of this program, dozens of SNe II have been 
monitored in the past ten years. We select the SNe II-P from the KAIT sample, since other core-collapse SN subtypes 
are not necessarily expected to obey the same correlations. However, the SN II-P subclass is not well defined. While 
histo rically SNe we re classified by their light-curve shapes, and SNe with a distinctive plateau were dubbed SNe II-P 
(e.g.. iBarbon et all 1979; Doggctt & Branch 1985), it is unclear how constant a plateau must be, how long that phase 
must last, and in which photometric bands it appears. In fact, those SNe which show hydrogen in their spectra but 
photometrically decline linearly (in magnitudes) were call ed SNe II-L, but not a single one of the best-studied objects 
(e.g., SN1979C) is considered a "normal" prototype (e.g. iPoznanski et al.ll200^ and references therein; Poznanski et 
al. 2009, in prep.). 

With the modern emphasis shifting to spectroscopic follow-up observations, the classification becomes, dishearten- 
ingly, even m uddier. It is o ften assumed by SN observers that only SNe II-P exhibit hydrogen lines having P-Cygni 
profiles (e.g.. ISchlegemi99"6[ ). but this claim has not been explored systematically. Lately, seve r al core-collapse SNe 
with extrenie luminosities have been discovered (e.g., lOfek et aLll2007t ISmith et al.l 120071 l2008t iQuimbv "eTaLl [20071 : 
iMiller et al.ll2008t iGezari et al.ll2008D . stretching the classification scheme beyond its current limits. 
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TABLE 1 
KAIT Sample 



lAU Name 


Heliocentric cz 


Flow cz^ 


cxpl 




{km s^"*") 


(kms^i) 


SN1999bg° 


1275 


1214(187) 


1259(5) 


SN1999em 


717 


819( 70)'' 


1476(4) 


SN1999gi 


592 


643(187) 


1520(4) 


SN2000bs 


8387 


8387(187) 


1650(6) 


SN2000dc= 


3117 


3113(187) 


1762(4) 


SN2000dj 


4629 


4625(187) 


1789(7) 


SN 2001x 


1480 


1329(187) 


1964(5) 


SN2001bq 


2623 


2238(187) 


2034(6) 


SN2001cm 


3412 


3413(187) 


2064(0) 


SN2001cy° 


4478 


4473(187) 


2086(6) 


SN2001do= 


3124 


3121(187) 


2134(2) 


SN2002gd"' 


2674 


2418(187) 


2551(2) 


SN2002hh 


48 


297(187) 


2576(2) 


SN2003hl 


2472 


2105(187) 


2868(4) 


SN2003iq 


2472 


2105(187) 


2921(1) 


SN2004du= 


5025 


5020(187) 


3223(4) 


SN2004et 


48 


297(187) 


3272(0) 


SN2005ay 


809 


816(187) 


3453(4) 


SN2005CS 


463 


505(187) 


3550(l)s 



(km s^-"-) 



nil 
(mag) 



y - /d 

(mag) 



Discovery Reference 



4657(439) 
3382(171) 
3697(200) 
4115(276) 
4545(232) 
4784(250) 
3800(203) 
4023(238) 
4065(189) 
3987(218) 
3927(311) 
2774(315) 
3692(435) 
3561(293) 
4043(179) 
4855(302) 
3980(194) 
3432(224) 
2102(211) 



15.31(0.02) 
13.23(0.03) 
14.00(0.02) 
17.67(0.09) 
16.05(0.06) 
17.00(0.09) 
14.74(0.02) 
15.55(0.03) 
16.30(0.02) 
16.13(0.04) 
15.75(0.03) 
16.78(0.03) 
13.68(0.05) 
15.60(0.02) 
15.27(0.03) 
16.60(0.04) 
11.39(0.01) 
14.67(0.02) 
13.98(0.02) 



0.74(0.04) 
0.58(0.03) 
0.86(0.03) 
0.51(0.10) 
0.97(0.06) 
0.59(0.11) 
0.57(0.03) 
0.51(0.08) 
0.68(0.03) 
0.63(0.05) 
0.97(0.05) 
0.74(0.05) 
2.67(0.05) 
0.98(0.04) 
0.66(0.03) 
0.53(0.07) 
0.52(0.01) 
0.63(0.06) 
0.69(0.03) 




Nakano et al. (2001) 
Jiang & Qiu (2001) 
Ganeshalingam et al. (20W) 
Modiaz & Li (2001)" 
Klotz et ak (2002) 
Li (2002) 
Li (2003) 
Llapasset (2003) 
Singer & Li (_2004) 
Zwittcr et al. f2004 ') 
Rich (20057 " 
Kloehr et al. (20051 



^ Corrected velocity of recession using the local Hubble-flow model o f ITonrv et al unless otherwise noted. 

^ Explosion date, assumed to be the midpoint epoch (given as Julian Day — 2,450,000) between the time of discovery and the 

last nondctcction (the uncertainty being half this value), unless otherwise noted. 

^ Fit velocity of the Fc II line at day 50, on the plateau. 

d Measured at day 50, on the plateau. 

^ Not used in the final sample, see ^4.11 

^ From the Cepheid distance of lLeonard et ahl ( 120031) and the Hubble law. 
S From lPastorello et ahl (1200611 . 

For the purpose of this work, we define the SNII-P subclass as follows: (a) prominent hydrogen features in the 
spectra, (b) no narrow emission lines indicative of interaction with the circumburst medium, and (c) a prominent 
"plateau" phase in the J-band light curve. 

The light-curve shape criterion is intentionally defined here rather loosely, so as not to preselect the best objects 
(beyond unavoidable observational biases) and is revisited in detail in t j4.1l Several dozen SNe in the KAIT sample pass 
these criteria, but for 19 of them we have the data necessary for empirical calibration and distance measurement; they 
constitute the main sample analyzed here (Table [1]). The SNe II-P that were not selected for this work were usually 
missing early-phase photometry to constrain the explosion date, or lacked spectroscopic coverage during the plateau, 
as required for measuring velocities. Two of our 19 ch osen SNe (SN 1999em and SN 1999gi) are already present in 
the H PN sample; they are very well-studied objects fe.g.. lHamuv et aI1l2001l : [Leonard et al.ll2002al lbl: iDessart fc HillieH 
l2006f ). We have used them to ascertain that there are no systematic offsets between our magnitude and velocity 
measurements and those of previous authors. Thus, our sample consists of 17 new SNe II-P. 

The optical CCD images of the SNe were reduced as follows. Flat-fielding and bias subtraction were performed 
automatically at the telescope using calibration frames appropriate for the science images. Galaxy subtraction and 
differential photometry were done using the KAIT pipeline (Ganeshalingam et al. 2009, in prep.). Galaxy template 
images were obtained with KAIT at least a few months after the SN had faded beyond detection. To ensure high-quality 
subtractions, templates were acquired on photometric nights with seeing < 3.0" . Two independent rou tines were used 
to perform galaxy subtraction. The first method is based on the ISIS package (lAlard fc Luptonlll998[ ) as modified by 
Brian P. Schmidt for the Hi gh-z Supernova Search Team (jSchmidt et al.|[l998D . The second method is based on the 
IRApi'' task PSFMATCH (IPhilli ps fc David I1995D . We performed differential point-spread function (PSF) fitting 
photometry to the results of both subtraction methods using the DAOPHOT package in IRAF to measure the SN 
fiux relative to local standards in the field. The results of the two subtraction methods were averaged. 

Calibrations were obtained on photometric nights using both KAIT and the 1 m Nickel telescope at Lick Observatory. 
iLandoltl (|1992D standards were observed at a variety of airmasses on each night to derive a photometric solution that 
we then applied to local standards for each of the SN fields. Instrumental magnitudes were transformed to standard 
Johnson-Cousins magnitudes using color terms derived from the photometric solution of many Landolt standard-star 
calibrations (see iModiaz et al ll200ll and Li et al.ll2001bl for more details). The uncertainty in our subtraction and 
photometry pipeline is estimated by injecting artificial stars with the same magnitude and PSF as the SN into the 
original KAIT images and recovering them. The final uncertainty is taken to be the scatter in recovering 20 artificial 
stars added in quadrature w ith the calibration error. We correct the magnitudes for Galactic extinction using the 
maps of Schlcg el et all (|1998l ). 

IRAF; the Image Reduction and Analysis Facility is distributed by the National Optical Astronomy Observatories, which is operated by 
the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation 
(NSF). 
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Fig. 1. — A representative spectrum (the closest to day 50) for each of the 19 SNe in our sample. The analysis herein makes use of about 
75 spectra of these 19 events. 

We obtained optical spectra with the Kast double spectrograph (jMiller fc Stonelll993D moun ted on the Lick Obser 



Oke et aLlll995f) m ounted on the 10 m 
Faber et al.ll2003f) on the 10 m Keck II 



vatory 3 m Shane telescope, the Low Resolution Imaging Spectrometer (LRIS; 
Keck I telescope, and the Deep Imaging Multi-Object Spectrograph (DEIMOS; 
telesco pe. The position angle of the slit was generally aligned along the parallactic angle to reduce differential light 
losses (iFilippenk o 1982). 

All spectra were reduced using standard techniques (e.g.. lFolev et al.l l2003'). Routine CCD processing and spectrum 
extraction fo r the Kast an d LRIS data were completed with IRAF. These data were extracted with the optimal 
algorithm of iHornd (Il986l). CCD processin g of the DEIMOS spectra was performed with a modified version of 
the DEEP pipeline fe.g.. IWeiner et 311120051 ). This produced rectified, sky-subtracted two-dimensional spectra from 
which onc-dimensional spectra were then extracted optimally (Hornc 1986). The wavelength scale was derived from 
low-order polynomial fits to calibration-lamp spectra. Small wavelength shifts were then applied to the data after 
cross-correlating a template night-sky spectrum to the sky spectrum extracted near the SN position on the slit. Using 
custom routines, we fit spectrophotometric standard-star spectra to the data in order to flux calibrate the spectra and 
to remove telluric lines llWade fc Horndri988t iMatheson et al.l[2000l ) . Information regarding both our photometric and 
spectroscopic data (such as observing conditions, instrument, data reducer, etc.) was obtained from our SN database 
(SNDB). The SNDB uses the popular open-source software stack known as LAMP: the Linux operating system, the 
Apache webserver, the MySQL relational database management system, and the PHP server-side scripting language 
(Silverman et al. 2009, in prep.). Figure [T] shows one spectrum of each of the 19 chosen SNe. 

Figure presents the BVRI light curves of our 19 SNe, and shows the epochs at which spectra were obtained. In 
total, there are about 1500 photometric points and 75 spectra used in this study. A thorough analysis of these data 
is beyond the scope of this paper, and will be discussed by Poznanski et al. (2009, in prep.). However, we note the 
great diversity of plateau shapes, a subject to which we return in H4.ll 

From the light curves we determine, using linear interpolation where necessary, the apparent magnitude in the / band 
and the V — I color at 50 d past explosion. We set the explosion date conservatively as the midpoint between the last 
nondetection of the SN and its discovery date (and the uncertainty as half this interval), except for SN 2005cs where we 
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Fig. 2. — Light curves of the 19 SNe II-P from KAIT: B (blue dots), V (green upside-down triangles), R (red triangles), and I (black 
squares). Vertical line segments at the bottom of each panel mark the times for which we have spectra. Note the abundance of spectra, 
and the diversity in the SN light-curve shapes. 



use the constraint determined by iPastorello et all (2006). For many objects better constraints on the explosion date 
could be set by using additional techniques (e.g., fitting templates to the photometry and/or spectroscopy); however, 
the explosion-date uncertainties have little impact on our results, since neither the color nor the magnitude (or Fe II 
velocities, see below) vary significantly during the plateau phase. We have well-sampled light curves, and as can be 
seen in Figure [1] we have spectra near day 50 for most of our SNe. 

Since the KAIT filt er system differs somewhat from the Bessell curves used by N06, we also S'-correct 
(iStritzinger et al" 120021) the photometry to the standard Bessell BVRI curves taken from the database of 
iMoro fc Munaril "! 2000( 1. We do the correction by warping a standard day 50 template spectrum to the KAIT photom- 
etry, and measuring synthetic magnitudes using the Bessell filter response. This has no noticeable influence on our 
results, since the corrections are all smaller than 2% in the measured values. 

From the spectra we measure the apparent Fe II velocity, defined here as the minimum of the Fe II line with respect 
to 5169 A. Unfortunately, there is no "industry standard" framework to determine the velocities of such broad, often 
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Fig. 3. — Left: Example application of our SNID-based Fe II velocity measurement to a typical spectrum. A day 90 spectrum of 
SN 2001X is cross-correlated with a set of templates using SNID, and a velocity is derived from each template (black dots). The best-fitting 
templates are selected based on SNID rlap values (blue circles), and sigma clipped to reject outliers (red crosses). A weighted mean is 
calculated from the selected velocities. Right: For a given object (SN 2001X in this example), each velocity (black circles) is propagated to 
day 50 (red circles) using Equation (2) of N06 (dashed line; our Equation jlj). A weighted mean supplies the final Fe II velocity for that 
SN, and its uncertainty. An additional uncertainty of 150kms~^ is added in quadrature to every SN, to account for peculiar motion. 

noisy, and sometimes asymmetric lines; the different approaches we tried reveal systematic differences in the derived 
values and hence different calibration parameters. In an effort to develop a method that is simple, robu st, and readily 
available to the community, we decided to use the supernova identification code (SNID; IB londin fc Tonrvl [20071. 
originally designed for the automatic spectral classification of SNc. It uses the lTonrv k David (l979l) cross-correlation 
algorithm, and is flexible enough to be used for our purpose. 

Figure [3] shows an example application of our method. We use the high-quality SNII-P templates distributed with 
SNID, for which the Fe II absorption-minimum wavelengths can be measured in a straightforward way (e.g., by fitting 
a Gaussian to the line) . We then cross-correlate every spectrum with the templates (using only the relevant wavelength 
range, 4500-5500 A), choose only those templates that produce a good fit (those with SNID rlap values greater than 
the median for the group), and calculate a sigma-clipped (Sti) weighted mean (left panel of Fig. [3]). We convert 
the wavelength offsets to velocity using the relativistic Dopplcr formula. Each velocity (and its uncertainty) is then 
propagated to day 50 on the plateau using Equation (2) of NOG, 



WFcIl(50d) =i;FeIl(t)(t/50d) 



0.464±0.017 



(1) 



Since most SNe have more than one available spectrum, we calculate for each SN the weighted mean of all derived 
Fe II velocities (right panel of Fig. [3]). We add to the velocity uncertai nty of every SN a valu e of CTpec =150 km s^^, in 
quadrature, to account for unknown peculiar velocities within the hosts (jSofue &: Rubinl[2001[ ). We have also confirmed 
that Equation ([T]) is consistent with our spectra, and we find no substantial residual from the NOG result, within its 
uncertainty. 

As for the HPN sample, the ap parent helioc e ntric redshifts of the SN hosts are derived from NED^-'^ and corrected 
for local flows using the model of iTonry et al.1 (|2000l ). The properties of SNe in our sample, as measured above, are 
included in Table [TJ 



2.3. Sample Comparison 

To understand whether HPN and KAIT sample a similar population, we compare the relevant observables in both 
samples. Figure S] shows a comparison of the distribution of uncorrected absolute /-band magnitudes (top panel), 
V — I colors (middle), and Fe II velocities (bottom), all at day 50 on the plateau, as measured above. These are 
the quantities used for standardization and distance measurement. Surprisingly, the luminosity distributions appear 
consistent with each other, while the KAIT sample features redder colors indicating a higher average extinction, and 
lower velocities that point to lower luminosities. 

Using Kolmogorov-Smirnov (KS) statistics, we find that the luminosity distributions match reasonably well (Pk-s 
values of O.GG), whereas the colors and velocities are less consistent (Pk-s values of 0.13 and 0.06, respectively). Those 
results remain similar if one considers only z < 0.1 SNe. However, the high- velocity objects in the HPN sample all have 
substantial uncertainties that are not taken into account by KS statistics. For example, most of the highest-velocity 
SNe in the HPN samples have velocity uncertainties of 2000 km s""'^. Nevertheless, differences in sample properties 
should be expected, due to the era in which they were compiled. KAIT is expected to be sensitive to fainter SNe, 
either intrinsically (hence with lower average velocities) or due to dust obscuration (hence with red colors). One should 

The NASA/IPAC Extragalactic Database is operated by the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 
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Fig. 4. — Comparison of the HPN and KAIT samples, showing observables relevant for distance measurements. Uncalibrated absolute 
magnitude (i.e., without applying dust or Fe II velocity correction; top), V — I (middle), and Fe II velocity (bottom). KS statistics indicate 
that the samples can be assumed to be drawn from the same underlying distribution at high significance for the absolute magnitude, but 
quite low significance for the color and Fe II velocities. 

note that while the KAIT sample may be less biased than the HPN sample, the luminosity function is still skewed 
significantly toward bright objects, due to follow-up criteria, as the median absolute magnitude is 0.5—1 mag brighter 
than the observed distribution of all KAIT SNcII-P (i.e., including those SNe that were not scheduled for follow-up 
observations; Li et al. 2009, in prep.). While the importance of this bias should be explored with future more balanced 
samples, we note that at high redshift the same selection naturally occurs. The difference between the samples in 
terms of the observed properties has a visible effect on the Hubble diagram we construct, as discussed in i i4.2l 

With an absolute magnitude distribution spanning more than 5 mag and a diverse set of light-curve shapes (Fig. [2]), 
the standardizability of SNe II-P would appear to be a significant challenge. However, as we demonstrate below, the 
resultant Hubble diagram remains tight, and the velocity vs. luminosity correlation persists, despite the larger and 
less homogeneous sample. 

Using the 23 SNe in the HPN sample, and the 19 in the KAIT sample, we construct our list of 40 SNe, after 
merging the two shared SNe, SN I999em and SN I999gi. For these two SNe we use the KAIT measurements, due to 
the somewhat smaller uncertainties (but consistent values) we obtain for the observables. Table [2] lists the full sample 
used in the subsequent analysis. 

3. FITTING METHOD 

NOG assume a correlation between luminosity and Fe II velocity that is only modified by a ICardelli et al.l (jI989l ) 
extinction law, with Ry = 3.1. However, there is no a priori reason to assume that all of the SNe in the sample are 
subject to the same extinction law, with the average curvature found in the Milky Way. For that reason we modify the 
color-dependent term in Equation (I) of NOG to be a free parameter. This term, which has a value of 1.36 in NOG^^, 
we refer to as i?/. Given the color excess EiV — I) = {V — I) — {V — /)o, Ri defines t he amount of extinct ion in the / 
band, such that Ai — RjE{V — I). The value of Rj can be directly derived from the ICardelli et all (|I989f) extinction 
curve as a function of Ry. We therefore obtain a slightly modified Equation (I) of NOG: 

Mi = Mi, - alogio(vFeii/5000) + Ri [{V - I) - {V - I)o] • (2) 

In addition, we use "Hubble-constant free" absolute magnitudes and luminosity distances Mig = Mi,— 5 \ogiQ{Ho) + 
25 and Vl = HqDl (|Goobar fc Perlmuttedfl995l : iSullivan et al.ll200Gn . since we use no anchor of "absolute" distance, 
and measure only the shape of the Hubble diagram (relative distances). We thus find 

A^/„ - a logio(z;Fo 11/5000) + i?/[(y-/) -0.53] -mj = -5log,o{'DL{z\n,n,^A)), (3) 

where we use for the ridge-line color the value from NOG, {V — I)o — 0.53 mag. As mentioned by NOG, the precise 
value for this color is degenerate with Mig, and does not infiuence our other results. As seen in the middle panel of 
fig. [4l perhaps due to the ubiquity of dust around SNe II-P, or as a consequence of an intrinsic dispersion, it is difficult 
to find a precise value for the unreddened color of those SNe. 

Since all of the terms in Equation ^ have errors, including significant redshift uncertainties for the nearby SNe, we 
maximize the marginal likelihood C by minimizing the multivariate cost function 



Note the sign misprint in N06 for that term. 
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TABLE 2 

Final Sample 



Name'^ 


czj z° 






y -r 


Best-fit \ji Source 




(kms"-'-) 


(kms^-"-) 


(mag) 


(mag) 


(mag) 



SN2002hh 

SN2004et 

SN2005CS 

SN1999gi 

SN1999br 

SN2005ay 

SN1999em 

SN1991G 

SN1992ba 

SN1989L 

SN1986I 

SN1990E 

SN2001X 

SN1990K 

SN2003hl 

SN2003iq 

SN2001bq 

SN1999ca 

SN2001cm 

SN1991al 

SN2000dj 

SN1992af 

SN2004dh 

SN1999cr 

SN1999eg 

SN1993A 

SN2000bs 

SN1993S 

SN1992am 

SN04D4fu 

SN03D4CW 

SN04Dlpj 

SN04Dlln 

SN03D3ce 

Culled SNe 

SN1999bg 
SN2002gd 
SN2000dc 
SN2001do 
SN2001cy 
SN2004du 



297(187) 
297(187) 
505(187) 
643(187) 
757(187) 
816(187) 
819( 70)d 
1029(187) 
1064(187) 
1189(187) 
1190(187) 
1273(187) 
1329(187) 
1623(187) 
2105(187) 
2105(187) 
2238(187) 
2772(187) 
0.011(0.001) 
0.013(0.001) 
0.015(0.001) 
0.016(0.001) 
0.017(0.001) 
0.019(0.001) 
0.019(0.001) 
0.027(0.001) 
0.028(0.001) 
0.029(0.001) 
0.042(0.001) 
0.119(0.001) 
0.138(0.001) 
0.139(0.001) 
0.185(0.001) 
0.257(0.001) 



1214(187) 
2418(187) 
0.010(0.001) 
0.010(0.001) 
0.015(0.001) 
0.017(0.001) 



3692( 435) 
3980( 194) 
2102( 211) 
3697( 200) 
1545( 300) 
3432( 224) 
3382( 171) 
3347( 500) 
3523( 300) 
3529( 300) 
3623( 300) 
5324( 300) 
3800( 203) 
6142(2000) 
3561( 293) 
4043( 179) 
4023( 238) 
5353(2000) 
4065( 189) 
7330(2000) 
4784( 250) 
5322(2000) 
4990( 300) 
4389( 300) 
4012( 300) 
4290( 300) 
4115( 276) 
4569( 300) 
7868( 300) 
3861( 150) 
3067( 660) 
4981( 214) 
3593( 159) 
5762( 522) 



4657( 439) 
2774( 315) 
4545( 232) 
3927( 311) 
3987( 218) 
4855( 302) 



13.68(0.05) 
11.39(0.01) 
13.98(0.02) 
14.00(0.02) 
16.67(0.05) 
14.67(0.02) 
13.23(0.03) 
15.01(0.09) 
14.65(0.05) 
14.47(0.05) 
13.98(0.09) 
14.51(0.20) 
14.74(0.02) 
13.87(0.05) 
15.60(0.02) 
15.27(0.03) 
15.55(0.03) 
15.56(0.05) 
16.30(0.02) 
16.06(0.05) 
17.00(0.09) 
16.46(0.20) 
17.60(0.15)== 
17.44(0.05) 
17.72(0.05) 
18.56(0.05) 
17.67(0.17) 
18.22(0.05) 
17.90(0.05) 
21.99(0.04) 
22.33(0.09) 
21.99(0.04) 
22.79(0.05) 
23.45(0.50) 



15.31(0.02) 
16.78(0.03) 
16.05(0.06) 
15.75(0.03) 
16.13(0.04) 
16.60(0.04) 



2.67(0.05) 
0.52(0.01) 
0.69(0.03) 
0.86(0.03) 
0.84(0.07) 
0.63(0.06) 
0.58(0.03) 
0.45(0.11) 
0.59(0.07) 
0.88(0.07) 
0.45(0.22) 
1.31(0.28) 
0.57(0.03) 
0.58(0.21) 
0.98(0.04) 
0.66(0.03) 
0.51(0.08) 
0.73(0.07) 
0.68(0.03) 
0.39(0.07) 
0.59(0.11) 
0.43(0.28) 
0.70(0.21)== 
0.56(0.07) 
0.55(0.07) 
0.51(0.07) 
0.51(0.20) 
0.69(0.07) 
0.38(0.07) 
0.59(0.05) 
0.90(0.10) 
0.66(0.06) 
0.70(0.07) 
0.39(0.50) 



0.74(0.04) 
0.74(0.05) 
0.97(0.06) 
0.97(0.05) 
0.63(0.05) 
0.53(0.07) 



28.89(0.23) 
28.36(0.09) 
29.61(0.21) 
30.57(0.13) 
31.60(0.43) 
31.27(0.13) 
29.84(0.11) 
31.70(0.36) 
31.33(0.18) 
30.94(0.19) 
30.82(0.29) 
31.44(0.37) 
31.59(0.11) 
31.62(0.75) 
32.01(0.14) 
32.16(0.10) 
32.55(0.13) 
32.93(0.89) 
33.18(0.10) 
34.28(0.54) 
34.26(0.18) 
34.05(0.87) 
34.86(0.30) 
34.56(0.18) 
34.68(0.16) 
35.68(0.17) 
34.71(0.31) 
35.32(0.14) 
36.27(0.11) 
38.85(0.11) 
38.52(0.47) 
39.28(0.11) 
39.43(0.12) 
41.22(0.79) 



32.41(0.18) 
32.90(0.21) 
32.93(0.14) 
32.35(0.15) 
33.01(0.12) 
33.94(0.13) 



KAIT 

KAIT 

KAIT 

KAIT 

N06 

KAIT 

KAIT 

N06 

N06 

N06 

N06 

N06 

KAIT 

N06 

KAIT 

KAIT 

KAIT 

N06 

KAIT 

N06 

KAIT 

N06 

N06 

N06 

N06 

N06 

KAIT 

N06 

N06 

N06 

N06 

N06 

N06 

N06 



KAIT 
KAIT 
KAIT 
KAIT 
KAIT 
KAIT 



^ Sorted by increasing redshift. 

^ From local Hubblc-flow model of ITonrv et al.l | |2000D (unless otherwise noted). For cz > 3000 
kms~"^ the redshift is given. 
== At day 50, on the plateau. 

^ From the Cepheid distance of lLeonard et all ll2003ft and the Hubble law. 

^ Note that these numbers differ from those in N06, following our re-reduction of the data. 



-2 log(Z:) ^ ^{ [-^/o - «logio(i'Fe 11/5000) + Ri \{V -I)- 0.53] - m, + 51ogio(PL(^))]' 

"'^logio(t;Fcll/5000) + ^i^{V-l) + ^m/ + ^5 log^ (Pi, (z)) + ^sys (4) 

+ log(Q;(Tfj,g^^(^p^jj/50Q(,) + Rjafy_j~i + a^^^ + crliog^^^VLiz)) + ^sys) }i 

where the summation is over the SNe in the sample (e.g., lKellvl[2"007[ ). The second, logarithmic term — which comes 
from the normalization of the likelihood function — helps mitigate the tendency of the first term to overfit the data 
by preferring large values of a, i?/, and agys- We note, however, that this term does not significantly change the 
resulting best-fit parameters. We include in the denominator (and the logarithmic term) of Equation ^ a systematic 
uncertainty agys in order to include (and measure) the intrinsic scatter in the correlation (i.e., the contribution to the 
dispersion not accounted for by the measurement errors). 

The parameters to be fit are therefore A4jg (although for simplicity we will keep citing values for A//„, assuming 
Hq = 70kms~^ Mpc~^), a, Rj (that we translate to i?y), cTsys, and the cosmological parameters. For the present low-z 
data, since we do not have any leverage on the cosmological parameters, we assume a standard cosmology {flm — 0.3, 
Q\ = 0.7), and only fit for crgys and the correlation coefficients. 
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Fig. 5. — Offset from the Hubble law fit to the full sample before culling, as a function of the /-band decline rate. As can be clearly seen, 
SN 2001cy and SN 2001do are amongst the strongest decliners and are our most significant outliers, followed by SN 1999bg which declines 
as well. The magenta dashed line marks our cutoff, where we select only SNe that have decline rates smaller than (or consistent to Icr 
with) zero. 

4. RESULTS 

Using the complete saniple of 40 SNe, we find that the best-fit values after marginalization are a = 4.6 ± 0.7, 
Rv = 1.5 ±0.5 {Ri = OJtti), Mi,, = -17.43 ±0.10 mag, and (Tsys = 0.38 mag. The systematic scatter is equivalent to 
an 18% uncertainty in distance, largely dominating the error budget, being more than twice the median measurement 
uncertainty. However, about a third of that scatter is due to two SNe in the KAIT sample, SN2001cy and SN2001do, 
and to some lesser extent SN 1999bg. Those SNe are offset from the best-fit solution by 6.7(7, 4.1(7, and 3.6ct, respec- 
tively. An analysis excluding those objects gives a tight correlation with Osys — 0.22 mag, which corresponds to 10% 
in distance. 



4.1. Sample Culling for Cosmology 

A possible glimpse at the source of the discrepancy lies in the photometric behavior of the outlying SNe. This is 
portrayed in Figure \5[ where we show, for the KAIT SNe, the distance from the Hubble law (as derived from the 
complete sample) in a units, as a function of the decline rate in the / band. We define here the decline rate as the 
slope of a linear fit to the /-band light curve starting 10 d after explosion (after the rise time), and ending at day 50, 
in units of mag (50 d)~^. The errors are estimated using a Monte-Carlo simulation propagating the uncertainties in 
the photometry and in the explosion dates. 

Figure [5] shows that SN2001cy and SN2001do have substantial decline rates, about 0.5 mag (50 d)^^, while most 
objects barely decline, or even have rising light curves. The decline vs. offset correlation is not strong; consequently, 
adding a decline parameter to the fitting only marginally improves the scatter to agys — 0.33 mag. However, we can 
now define quantitatively a selective subsample of SNe II-P for the purpose of distance measurement in this paper, 
and reject objects in a controlled manner. For the purpose of this work, we will keep objects that have decline rates 
smaller than (or consistent to la with) zero.^^ 

As can be seen in Figure[5l this criterion rejects all three outliers (SNe2001cy, 2001do, and 1999bg) by design, but 
also SNe2000dc, 2002gd, and 2004du that are relatively well fit. This choice leaves 34 SNe in the combined sample 
that we analyze in the following sections. We emphasize that a similar cut has been applied by N06 prior to the 

Note that we are defining the specific subsample of SNe II-P in this way only for the purpose of having a clear and clean set for 
cosmological purposes. We do not mean to imply that every SN II having a decline rate larger than zero is necessarily a SN II-L or some 
other "not-II-P" SN. The real distinction between SNe II-P, SNe II-L, and other potential SNII subtypes will be considered in more detail 
by Poznanski et al. (2009, in prep.). 
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Rv Rv 

Fig. 6. — Best-fit parameter {a, Mj^^, Rv) lo", 2(t, and 3o" contours at left; marginaUzed posterior probabilities at right, for the culled 
sample of 34 SNe. 

costly follow-up observations of the moderate-z SNe in their sample (P. Nugent, 2008, private comm.), and has most 
probably been applied to the HP02 sample as well. The precise criteria have not been detailed previously, and could 
significantly bias the best-fit parameters and dispersion. Our choice, while somewhat arbitrary, does not affect our 
results below, and any choice that rejects the two major outliers will effectively be equivalent. 

4.2. Best-Fit Parameters and Hubble Diagram 

Using the culled sample of 34 SNe, we find that the best-fit values after marginalization are a — 4.4 ± 0.6, Ry — 
1.7±0.5 (i?/ = 0.8±0.3), M/„ = —17.39 ±0.08 mag, and CTsys = 0.22 mag (10% in distance), consistent with the values 
for the complete sample. The likelihood contours and the marginalized posteriors are plotted in Figure [S] As can be 
seen, the coefficients are weakly covariant. 

The value of the best-fit aver age Ry, while small, is consistent with results for many SNe la, which tend to suffer an 
extinction that is very s elective ()Ehas-Ros a et al." 2006HKris ciunas et al.n 2007l : lElias-Rosa et al]|2008l : [Nobili fc Goobarj 
l2008HWang et al.ll20'08f ). As shown by|Wang (2003) and|Goobar (2008), this result could be explained with a fairly 
simple model where "normal," yet circumstellar, dust causes heavily color-dependent scattering. While this model is 
discussed in the context of SNe la, core-collapse SNe are expected (and observed) to be enshrouded by more circum- 
stellar dust, which could explain these low Rv values. We note that a small Rv is still the best fit if one removes the 
highly extinguished SN 2002hh from the analysis. 

Using Equation ([3]) we build the Hubble diagram in Figure [71 The scatter around the standard cosmology trend is 
quite small, and is accounted for by the systematic and measurement uncertainties. The systematic scatter can be 
contained in a reasonable velocity dispersion of cFv-paii — 135 kms^^. This is equivalent to, or better than, the results 
obtained by previous authors, and is encouraging considering the size of the sample and its initial inhomogeneity. 

4.3. Robustness of the Correlation 

In order to test the robustness of the correlation, we explore the importance of the different ingredients in the fit. If 
one applies no correction at all (i.e., assume the SNelTP are perfect standard candles), the correlation crumbles, and 
gives (Tsys ~ 0.9 mag. The same occurs if one does not apply the velocity correction, but only the color correction. 
A simple test of the importance of that term is obtained by shuffling the velocity measurements among the objects, 
assigning to each SN a random value chosen from the sample distribution. The result, on average, is a very weak 
correlation, w ith a^ys ~ 0.6 mag, and unrealistic best values for the coefficients (a < and Rv > 15 in most cases; see 
lDrainel[2003l for a review on measured values and theoretical limits on Rv)- 
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Fig. 7. — Hubble diagram (top), and residuals from standard concordance cosmology (bottom), for the HPN sample (black triangles) 
and the culled KAIT sample (orange filled circles). "Standard" (line) and Einstein-deSitter (dash-dot) cosmologies are added to guide the 
eye. The bottom-right panel shows the similarity in the distribution of residuals for both samples. 

However, if we use solely the Fe II correction (which is numerically equivalent to setting i?/ = or Ry = 0.78), we 
get a scatter only slightly greater than for the best solution (csys ~ 0.26 mag), and mildly different values for a and 
Mjg. This result indicates that for our sample of SNe, the dust correction is not strongly required by the data. 

Since we do expect SN magnitudes to be affected by dust, there are a few possible explanations we consider. First, it 
could be that an intrinsic color- velocity correlation masks most of the contribution from dust. While theoretically one 
could expect the color (i.e., the temperature) to be correlated with the photospheric velocity at least to some extent, 
we find no indication in the sample for any such covariance, despite having at least a few SNe suffering negligible 
extinctions. 

An alternative explanation is that the sample is heavily biased toward dust-free objects. While this is se curely 
wrong for at least one object (SN2002hh, which suffers ~5 mag of extinction in the V band; iPozzo et al.l 120061 ). it is 
probably wrong for many of the other SNe as well. HG5 finds significant dust corrections for at least some of the SNe 
in his sample, and our sample should be less "hand-picked" as the objects included were discovered by modern-era 
CCD-equipped SN searches that can find more SNe II-P buried in their host galaxies. 

However, neglecting SN2002hh, most SNe in the sample are within ^ 0.3 mag of the sample's mean V — I color. 
Additionally, for any reasonable dust law, the color term in Equation ^ is at least 3-5 times less significant than the 
velocity term, so that even substantial differences in color will contribute relatively little to the dispersion in distance 
moduli. Consequently, a sample that does not have many heavily extinguished objects can be fit about as well when 
assuming no extinction at all. 

We have searched unsuccessfully for parameters that correlate with the residuals from the Hubble diagram and 
further reduce the necessary cTgys value, other than the slope of the plateau in the light curve that we have used to cull 
our sample in §4.11 The i3-band luminosity at day 50, for example, does not reduce the scatter by more than a few 
percent of its previous value. There are indications that rejection of red objects (those having {V ~ I) ~ {V — I)q > 0.2 
mag) may reduce the scatter, but our sample is still too small for us to make a robust statement regarding this. 

We also determine the best-fit parameters for different subsamples of the data. The HPN and KAIT samples when 
analyzed separately give best-fit parameters consistent to within la with those derived from the full set. When 
examining various cuts in redshift, we find that the lowest-redshift objects in the sample {z < 0.004) tend to pull the 
solutions to somewhat larger values of a (near 6 instead of 5), as previously noted by N06 and tentatively ascribed to a 
Malniquist bias. At z > 0.004 there seems to be a systematic difference between the two samples, with the KAIT SNe 
being mostly underneath the Hubble-law line. The number of SNe in this range is too small for a conclusive analysis, 
but we estimate that this is a reflection of the differences between the samples noted in ^2.31 KAIT finds intrinsically 
fainter and more extinguished SNe, for which the fit compensates by reducing their derived distances. 



4.4. Shared-Host SNe 

SN 2002hh and SN 2004et both occurred in the same host galaxy, NGC 6946, a neighbor of the Milky Way (< 10 Mpc 
away) . While SN 2004et shows no apparent uniqueness (except perhaps some small photometric jitter in the plateau 
phase of the li ght curve), SN 2GG2hh is highly redden ed by dust, with Ay ~ 5 mag and an infrared echo from a shell 
with ~10 M0 1 B arlow et al . 2005; P ozzo et al.l 1200(1 ) . Despite the extreme extinction, the distances of the two SNe, 
as derived using the best- fit parameters of M.2\ are consistent within about 2a (see Table [2|). 

SN 2003hl and SN 2GG3iq exploded within weeks of each other, in the same host galaxy, NGC 772 (see F ig. jSl 
left panel). In fact, SN2G03iq was discovered by an amateur SN observer while following 2003hl (jLlapassetl [20G3i) . 
With the exception of the first spectrum of SN 2003hl, our spectra were obtained by placing the slit on both objects 
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Fig. 8.— Le/t; Color composite image of NGC 772 showing SN2003hl and SN 2003iq. It consists of KAIT BVRI images, and a deeper 
frame from Deep-Sky (Nugent et al. 2009, in prep.) not showing the SNe, but enhancing galaxy features. Right: Distance difference 
between SNe 2003hl and 2003iq (solid line) and between SNe 2002hh and 2004et (dashed line), in a units, as a function of the assumed 
Rv- A value of Ry =3.1 (dashed grey line) is rejected at more than Aa. The gray area marks the Icr interval preferred by the full sample. 
These two pairs of SNe, as well as the full sample, favor a low-ijy dust model. 

simultaneously. In this lower-extinction case (compared with SN 2002hh in NGC 6946) , the agreement in distance is 
even better. 

However, the distances to both pa irs of SNe be st agree only for the favored model with low Ry- (A similar value 
for SN2002hh has been measured bv lPozzo et~anr2006 using a full light-curve comparison to SN 1999em.) As seen in 
the right-hand panel of Figure [H a value of Ry = 3.1 is rejected at a combined significance level higher than 4a. This 
result strongly supports the low Ry value favored by the full sample. 



/Rv = 1.7 ±0.5, 
/ best-fit 
using 

full sample 



4.5. Error Budget: Prospects for Future Samples 

While the scatter we find in §4.21 is small, it is tightly coupled to our error estimates. The value agys = 0.22 mag 
is in effect the scatter unaccounted for by other error terms. As can be seen in Table [21 due to the quality of our 
data, our uncertainties are generally smaller than those of previous authors. If one argues that those uncertainties are 
underestimated (though we find that they are not), increasing them by some amount will result in an even smaller 
value of (Tsys- 

As a consequence of our precise measurements, systematic errors dominate our Hubble diagram. This allows us to 
constrain the intrinsic dispersion of the correlation to be roughly 10% in distance. However, one should note that 
there are observational limitations that would increase that uncertainty for any foreseeable sample. Even for data 
with exquisite photometry, and high S/N spectroscopy, there still remains an uncertainty in the peculiar velocity of 
the exploding star, on the order of 150kms~^. This translates to about 0.05— 0.08 mag of uncertainty, i.e., 2—4% in 
distance. A reasonable precision to expect for a SN with a single spectrum is about 300kms~^, or ~0.12mag. 

In addition, any project that wishes to amass a suitable sample of SNell-P will need to address the critical issue of 
sample definition discussed in ij4.1l A simple strategy to overcome this would include high-cadence observations during 
a "detection phase," in order to have strong constraints on the explosion dates, followed by low-cadence monitoring of 
the candidates, in order to ascertain that they do not decline. 

The weak dependence onV — I color discussed in §4.31 has strong utility for extension to high-z SNe, where selection 
against faint and reddened SNe is stronger, and dust could have less impact. A Hubble diagram based only on the 
V band, without any color term, is equivalent to fixing i?/ = — 1 in Equation ([3]). Solving for the other parameters, 
we find a somewhat inferior scatter of (Tgys — 0.29 mag. However, obtaining rest-frame /-band magnitudes for z > 0.4 
SNe is observationally challenging, requiring near- infrared observations. The slightly poorer scatter in the V band 
can be overcome by statistic s. We do no t expe ct significant scatter (or evolution) in those photometric bands due 
to metallicity. Figure 13 of iBaron et al.l (|2003l ) shows that wildly different metallicities produce little variance in 
the intrinsic colors of SNell-P at wavelengths longer than ~5000A. Future low-z samples will also allow a better 
quantitative determination of the cutoff in decline rate that separates "good standardizable SNell-P" from other SNe. 

5. CONCLUSIONS 

Using a sample that is larger and more diverse than those studied previously, we have shown that SNell-P can be 
calibrated using simple correlations, yielding a tight Hubble diagram with a dispersion of ^10% in distance, not much 
worse than for SNe la, if one carefully constructs the sample of SNe. We have shown that the N06 standardization 
method, a derivative of the HP02 method, might be simplified even further because the correction for dust extinction, 
based on the current sample, has a low statistical impact on the scatter in the Hubble diagram. This implies that 
a distance measurement could, in principle, be obtained with a single spectrum of a SN during the plateau phase, 
combined with rest-frame V or J-band photometry at a similar time. However, additional (lower S/N) photometry is 
required during the plateau phase, in order to reject non-II-P SNe. 
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If dust correction is applied, the best-fit solution prefers a very steep dust law, with Rv < 2, as recently indicated for 
SNe la. This may have significant implications for SNIa cosmology, where systematic uncertainties currently dominate. 
It is further supported by two pairs of SNc that occurred in the same host galaxies, and whose distances agree only 
for low Rv values. 

The Hubble diagram in Figure[7]is dominated by nearby objects, where systemic velocities of the host galaxies govern 
the error bars. More exact values of the parameters, and more confidence in the method, will require the analysis of a 
sample of SNe in the Hubble flow. The most useful redshift range to calibrate the methods is z « 0.05-0.1, a distance 
that has been almost inaccessible until recently due to the requirement of a very large, and moderately deep, search. 
Pan-STARRS f Ka iser et all 12002 ?) and the Palomar Transient Factory (PTF; Rau et al. 2009, in prep.) wiU supply 
more SNe II-P than any follow-up program could realistically handle. If subsequent follow-up resources are allocated, 
those SNe will become the backbone of any cosmological use of SNc II-P, anchoring the correlation on more secure 
grounds. 
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